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Abstract 



We investigate de Haas-van Alphen (dHvA) oscillations in the mixed state 
of a type-II two-dimensional superconductor within a self-consistent Gor'kov 
perturbation scheme. Assuming that the order parameter forms a vortex 
lattice we can calculate the expansion coefficients exactly to any order. We 
have tested the results of the perturbation theory to fourth and eight or- 
der against an exact numerical solution of the corresponding Bogoliubov-de 
Gennes equations. The perturbation theory is found to describe the onset of 
superconductivity well close to the transition point H C 2. Contrary to earlier 
calculations by other authors we do not find that the perturbative scheme 
predicts any maximum of the dHvA-oscillations below H C 2- Instead we ob- 
tain a substantial damping of the magnetic oscillations in the mixed state as 
compared to the normal state. We have examined the effect of an oscillatory 
chemical potential due to particle conservation and the effect of a finite Zee- 
man splitting. Furthermore we have investigated the recently debated issue 
of a possibility of a sign change of the fundamental harmonic of the magnetic 
oscillations. Our theory is compared with experiment and we have found good 
agreement. 
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PACS numbers: 74.25 Ha, 74.60-w 



I. INTRODUCTION 

In recent years there has been a renewed interest in the interplay between external mag- 
netic fields and superconductivity in type-II superconductors. It is well known that de 
Haas- van Alphen (dHvA) oscillations are a useful tool for probing the Fermi surface in met- 
als in the normal state. For type-II superconductors the magnetic field is allowed to partially 
penetrate the sample in the mixed state. One would then expect magnetic oscillations in 
the mixed state to give information about the quasi-particle dispersion and the magnetic 
field dependence of the correlations in the ground state. Magnetic oscillations in the mixed 
state were observed for the first time in the layered superconductor 2H — NbSe2 over 20 
years ago.0 More recently dHvA oscillations were observed in the organic superconductor 
k — (ET) 2 Cu(NCS) 2 $ the A15 compounds V3S2! and Nb 3 Sn$ the borocarbide supercon- 
ductor YNi 2 B 2 C$ and the high temperature superconductors Y BaCu(M and BaKBi(M. 
These experiments have sparked a variety of theoretical investigations, not least in order to 
understand the interplay between oscillations in the quasi-particle spectra and the ground 
state condensation energy. The transition line H c2 between the normal state and the mixed 
state was shown to exhibit weak oscillations as a function of the magnetic field.!'! For high 
magnetic fields, clean samples, and very low temperatures H c2 has been predicted theoreti- 
cally to be a strongly oscillating function.0 The mixed state is characterized by the interplay 
between Landau level quantization due to the magnetic field, and Cooper pair formation 
characteristic of superconductivity. This calls for a theory that takes both effects into ac- 
count consistently. The theory developed by Maki0 and StephenH gives a simple picture of 
the vortex lattice acting as an extra scattering potential on quasi-particles thereby damp- 
ing the magnetic oscillations. The theory uses semiclassical approximations and, crucially, 
fails to impose the physical condition that the vortex lattice is the self-consistent mean 
field of the Cooper pairs. The problem simplifies when the electrons are confined to form 
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pairs within the same Landau level (diagonal approximation) and this case has been treated 
by several author sJllllfl Unfortunately the diagonal approximation ignores the fact that the 
typical excitation is a superposition of an electron and a hole in different Landau levels, 
but with similar energies. This effect is strongest when the chemical potential \i is either 

at a Landau level rif = ^ 1/2 = n (n integer) or exactly between two Landau levels 

rif = n + 1/2. We then have exact degeneracy between an electron state in a Landau level 
rif + m and a hole in the level rif — m, when rif = n, and between an electron in a level 
rif + m + 1/2 and a hole in a level rif — m — 1/2, when rif = n + 1/2, respectively. A ma- 
jor effect of the self-consistent pairing field is then to mix these two degenerate excitations 
strongly. Following the results of the diagonal approximation Dukan et aZ.IHI have focused on 
the consequences of a gapless portion of the quasiparticle spectrum. The calculation, which 
is appropriate for low lying excitations in 3 dimensions, is not applicable for two dimensional 
systems where the number of gapless points and their dispersion law vary strongly with the 
magnetic field and it does not take into account the oscillatory behaviour of the ground state 
energy as a function of the magnetic field. This oscillatory behaviour of the ground state 
energy has been considered by P. Miller and B. L. Gyorffyii in the A 3> k^T limit. Norman 
et a/.0 have studied the problem numerically and have linked the damping of the magnetic 
oscillations to the broadening of the Landau levels due to the gap. Recently^ there has 
been claims based partly on Gor'kov theory and partly on an assumed simplified form for 
the quasiparticle spectrum that below a certain field H inv < H c2 , the magnetic oscillations 
should exhibit a rapid 180° phase shift. 

In this paper we develop a new scheme for calculating the Gor'kov expansion terms 
treating the quantum effects of the magnetic field exactly. In addition we solve numerically 
the corresponding Bogoliubov-de Gennes (BdG) equations. Using the developed formalism 
we study the magnetic oscillations in the mixed state of a type II superconductor. We are 
working in two dimensions since many organic metals are known to show almost perfect 2D 
behaviour. Exploiting the symmetry of the magnetic translation group of the vortex lattice 
we have been able to calculate the expansion coefficients in the Gor'kov theory exactly to any 
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order making no restriction on the energy of the center-of-mass of the Cooper pairs. Self- 
consistency within this approach then transforms to the much simpler problem of minimising 
a polynomial of a finite number of variables. This allows us to develop an analytical theory 
for the thermodynamic potential and thus for the magnetic oscillations close to H C 2 which 
contains no approximations apart from the assumption of a small order parameter. This 
establishes a rigorous basis for our theory, compared with earlier attempts. It turns out to 
be crucial to determine the order parameter self-consistently since its oscillatory behaviour 
when the magnetic field varies is the cause of the damping of the dHvA oscillations. We find 
that the dHvA oscillations are damped in the mixed state as compared to the normal state, in 
agreement with what is observed experimentally. This is due to the fact that the contribution 
from the superconducting order parameter to the magnetic oscillations partly cancels the 
contribution from the normal grand potential. The superconducting order parameter itself is 
an oscillating function of the magnetic field, with local maxima occurring whenever we have 
a Landau level at the chemical potential since electrons can then form Cooper pairs without 
any cost in kinetic energy. This is the simple physical picture of the damping emerging from 
our formalism and it complements the interpretation given by P. Miller and B. L. GyorfljEl 
and Norman et a/.0 When many Landau levels participate in the pairing we have simplified 
the expressions for the expansion parameters. This makes it possible to give fairly simple 
analytical expressions for the rate of damping af the dHvA oscillations close to the transition 
line that may prove useful when fitting experimental data. 

A similar approach has been taken by Maniv et aZ.EH Using the semiclassical and various 
other approximations, they calculate the Gor'kov expansion coefficients for a 2D metal to 
fourth order in A(r) when the motion of the centers of mass of the Cooper pairs is restricted 
to the lowest Landau level. However, they obtain^ that the magnitude of the magnetic 
oscillations exhibit a maximum below H c2 . This is contradicted by our exact calculation of 
the expansion coefficients and also by our numerical solution to the BdG-equations. 

Recently it has been suggested that the degeneracy of the Landau levels should give 
rise to non-perturbative terms in the expansion of the grand potential thereby making the 
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traditional Gor'kov theory invalidEil We have tested our perturbative theory carefully against 
an exact numerical solution of the BdG-equations and we do not find any of the predicted 
non-perturbative effects. The theory based on the Gor'kov expansion agrees very well with 
the exact solution if we are not too far below H C 2. It is essentially a high temperature 
expansion in the sense that is an asymptotic series as long as the change in the quasiparticle 
levels as compared to the normal state is not larger than ~ 0(&&T).§1 

In a two-dimensional metal the chemical potential is an oscillatory function of the mag- 
netic field when the number of particles N is fixed. When higher harmonics are important 
(i.e. low temperatures and clean samples) the dHvA signal in the normal state for fixed N 
look qualitatively different from the case when the chemical potential is fixed. It is of interest 
to see what consequences this difference has for the magnetic oscillations in the mixed state. 
Examination of the dHvA oscillations in the mixed state in the two cases yields that the 
superconducting order for fixed number of particles reduces the oscillations in the chemical 
potential and that the dHvA oscillations are essentially the same in the two cases apart from 
a narrow region close to H c2 . Specifically, the rate of damping of the magnetic oscillations 
is the same when the number of particles is constant and when the chemical potential is 
constant. 

Since the contribution to the magnetic oscillations from the condensation energy is in 
antiphase with the contribution from the normal grand potential, it has been suggested^ 
that this will result in a sign change of the fundamental harmonic of the dHvA oscillations for 
H < H inv < H c2 . This would happen if the superconducting contribution were to overwhelm 
the contribution from the normal grand potential deep enough into the superconducting 
state. Based on an approximate evaluation of the Gor'kov expansion parameters one can 
calculate an expression for H inv B Using our expressions for the damping, we are able to 
predict that within the region of validity of the perturbative scheme this effect will not occur. 
Hence there is no theoretical reason, within perturbation theory, to expect inversion of the 
magnetic oscillations. This result agrees with the lack of experimental observation of such 
an effect. It also agrees with our exact numerical solutions to the BdG equations which show 
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a complete suppression of the magnetic oscillations deep enough into the mixed state.il 

Although there are currently experimental uncertainties about the value of H c2 in the 
organic superconductors, a comparison with experimental results for the quasi 2D supercon- 
ductor n — (ET)2Cu(NCS)2 yields good agreement between theory and experiment. 

The outline of our paper is as follows. Sec. |TJ sets up the formalism for describing the 
vortex state using both the Bogoliubov-de Gennes equations and perturbation theory. In Sec. 
I [J we compare the results of the perturbation theory with the exact numerical solution. The 
damping of the magnetic oscillations is discussed in Sec. [TV[ We give a physical interpretation 
of the damping. The effect of a finite Zeeman splitting term is discussed and the case of a 
conserved number of particles as opposed to a conserved chemical potential is considered. 
Using approximate expressions for the damping parameters we are able to give a simple 
analytical expression for the rate of damping of the dHvA oscillations close to H C 2. The spin 
dependence and the temperature dependence of the oscillations can then be extracted. We 
then examine the validity of the arguments leading to a sign change of the first harmonic 
of the dHvA oscillations. In Sec. |VJ we compare our analytical theory with experimental 
results. Finally we summarize our results in Sec. |V11| . 

II. ELECTRONS IN THE VORTEX STATE 

A. General representation and BdG- equations 

We consider a pure 2D electron gas in the x—y plane with a perpendicular magnetic field 
H along the z-axis. In the Landau gauge, A = (0, Hx, 0), the single-particle eigenstates can 
be chosen to be 

1 T _ kl 2 

WW = ^=e-^<M^^) (1) 



where 4>n{.x) = (2 N Nly/nl)' 1 ^ 2 HN(x)e~2 x2 with being a Hermite polynomial of order N, 
and I 2 = hc/eH is the magnetic length. The size of the system is L x x L y . Band structure 
effects are assumed to be adequately described employing an electron effective mass m* . 



The B-field is taken to be uniform within the sample thereby ignoring the partial screening 
by the supercurrents. This approximation holds for strong type-II superconductors (k ^> 1) 
such as the organics, where the penetration depth is much larger than the coherence length. 

In the mixed state of a conventional type II superconductor the order parameter forms 
a vortex lattice. It is therefore advantageous to use a basis set which incorporates this 
symmetry. We have chosen to use the following set of functions introduced by Norman et 

^k(r) = x ^T,^ a ^ t2/ ^N,- ky+ taM^ ( 2 ) 

V ^ t 

where k x £ [0, with Ak x = |^ and k y e [0,a x /l 2 [ with Ak y = j^- define the magnetic 
Brillouin zone (MBZ). The symmetry of the order parameter restricts the pairing to be 
between electrons with quantum numbers k and -k.El By adjusting a x we can obtain both 
a triangular (a x = Z(\/37r/2) 1//2 ) and a square vortex lattice (a x = l(n/2) 1 ^ 2 ). Throughout 
this paper we choose to work with the triangular lattice since we expect the free energy to 
be minimized by this symmetry (except possibly in the re-entrant regime) .0 We are using 
mean field BCS-theory with a smooth cutoff in the interaction around the Fermi surface, 
applicable for weak-coupling superconductors. The mean-field Hamiltonian is: 

H = H + Ht 

H = J dv^liv) ( (P ~| A)2 -^) Mr) 

H x = Y,l d * [A(r) W (iV) W (M)0^ k (r)^_ k (r)a^ kT a] v _ ki + ex.] (3) 

NM J 

k 

where the order parameter is defined as: 

A(r) = VJ2 w{N)w{M)<f) Nk {r)<f) M _ k (r) < a Nk1 a M _ H > (4) 

ATM 

k 

This differs from the conventional BCS-hamiltonian since we have introduced the weight 
function w(n). It is necessary to have a smooth cutoff in the pairing interaction since we 
otherwise would get non-physical effects arising from Landau levels abruptly entering or 
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leaving the pairing region. The weight function w(n) is chosen to be Gaussian i.e. w(N) oc 
e -(^N/Q-5hw D ) 2 w j iere UJd j s the pairing width and £jy = {N + l/2)hu c — /i. This approach 
was introduced by Norman et although they used a different weighting function. It 
should be noted that the above slightly unconventional definition of the order parameter 
is necessary. Otherwise the self-consistency condition is not equivalent to minimising the 
grand potential f2 with respect to A(r) (i.e. g^rp, =0). In the vortex lattice case the order 
parameter can be characterized by a finite number of parameters A.,0 



The Aj-'s are determined selfconsistently as explained in reference 17. Assuming not only 
translational but also six fold rotational symmetry of |A(r)| gives the restriction^ j = 
0, 6, 12, . . . where j < 2iV max . iV max is the highest Landau level participating in the pairing. 
Using the above transformation the corresponding BdG-equations split into a set of equations 
for each k and they can be solved numerically. Norman et a/.0 have carried out an extensive 
numerical investigation of the quasiparticle spectrum and the magnetic oscillations in the 
superconducting state. We have developed a similar numerical scheme to solve the BdG- 
equations. In this way we are able to check our analytical results against an exact numerical 
solution. 



Since we are interested in the region near H c2 where the order parameter is small, it is 
natural to consider the Gor'kov expansion of the grand potential. This can be done either 
through the equation of motion approach originally used or by using the grand partition 
function for the symmetry-broken self-consistent Hamiltonian: 





B. Perturbative expansion of the grand potential 
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(6) 



S 



where T>(ip*(r, T)ip a {r, r)) denotes functional integration over Grassman variables. We have 
defined ip a (r) = I] njk w(n)</> nik (r)a nikcr . The -^|A(r)| 2 term corrects for the double counting 
of the interaction energy in the Hartree-Fock approximation. Expanding the grand potential 
Q = — -a In 2 in powers of A(r) we obtain to eighth order 

tt S - tt N = tt 2 + fi 4 + ^6 + ^8 (7) 

where 

Q 2 = dr|A(r)| 2 -±J dr 1 dr 2 A(r 1 )A*(r 2 )K 2 (r 1 , r 2 ) 
^4 = Jdv 1 . . . dr 4 K 4 (r h r 2 , r 3 , r 4 )A(r 1 )A(r 2 )A*(r 3 )A*(r 4 ) 

Qe = -^Jdr 1 ... dr 6 tf 6 (n, . . . , r 6 ) A(n) A(r 2 ) A(r 3 )A*(r 4 ) A*(r 5 )A*(r 6 ) 
fi 8 = ^ / dn . . . dr 8 tf 8 (n, . . . , r 8 ) A(n)A(r 2 ) A(r 3 )A(r 4 ) A*(r 5 )A*(r 6 ) A*(r 7 )A*(r 8 ) (8) 

The kernels are given by 

K 2 (r 1 , r 2 ) = i 53 <^(r 2 , ri, -a>„)G°(r 2 , r x , u^) 

K 4 (ri, r 2 , r 3 , r 4 ) = i 53 G |( r 4, r i7 w^G^ra, ri, -w„)GJ(r 3 , r 2 , ^)G^(r 4 , r 2 , -u^) 

if 6 (ri, • • • , r 6 ) = -g 2 G K r 6 ; r i ; ^)G , °(r 5 , n, -w„)G°(r 5 , r 2 , w„) 
xG°(r 4 , r 2 , -^)GJ(r 4 , r 3 , u v )G°(r 6 , r 3 , -a>„) 
i^ 8 (ri, . . . , r 8 ) = p 53 ^K r 6j r i; w„)G^(r 7 , ri, -w„)GJ(r 7 , r 2 , ^)(5°(r 8 , r 2 , -u v ) 

xGJ(r 5 , r 3 , u v )G°(r 6 , r 3 , -ov)G°(r 8 , r 4 , ^)(5°(r 5) r 4 , -u>„) (9) 

and c<j„ = (2v + l)7rkf,T/h are the Matsubara frequencies. Maniv et a/.El have calculated the 
expansion up to fourth order in A(r) using essentially semiclassical approximations. They 
used a variational form of the order parameter which has no symmetry built in initially but 
restricts the electrons to condense in the lowest center-of-mass Landau level (Aj^ = 0). 
As will be shown below, this restricion introduces no serious error within in the region of 
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interest in the phase diagram. Since it is known0 that the triangular lattice is the minimal 
energy configuration (except for the re-entrant regime) we have exploited this symmetry to 
calculate these expansion terms exactly. Because we are using a smooth pairing cutoff in 
our Hamiltonian we have, instead of the Green's function for the normal state G° a (r^i v ii ^u), 
the following function in our kernels: 

riot \ 0nk(l"2)0^k( r l) 2/ \ /in\ 

G a {r 2 , ri, u v ) = ^ w (n) (10) 

where £ n(T = £ n + gm*a/2motkv c . The only difference from the Green's function for the 
normal state is that we have included the weight functions w(n) in the sum. Using the 
symmetry of the vortex lattice the integrals can be solved. We have to fourth order 



and 



V2lL X Ly j 

E«A 2 (ii) 



V 4 a 4 

^4 = gr4r I 74 2 w2 ( n i) w2 (^2)w 2 (n3)w 2 (n4)/(ni,n2,n3,n 4 ) 

° L x L y L m...n 4 



E R mn 4 R n 3 n 2 R nm a R n3n4^ni+n4-ji,n2+n3-j2 A AAA 
il J2 -°J3 -°i4 ^n 1 +n 2 -js,n :i +n i -j4. LA Ji LA h L ^j3 l ^j 
1...J4 

E 7n . . .34 A? 1 Aj- 2 A 3 - 3 Aj 4 (12) 



J1-J4 



where 



and 



s fe£ = ^E^ 1J2 ^f 4 E e***-*® [<f> jl+j2 ; (2/^,k> ; ,, ;; ,(2/^.,.). 

-i7r(/ii-h, 2 ) 



e-^-^^+j^j^hidx + r/,.)o ; ,. ;; ; (2// 2 r/,. + a x )j (14) 

The coeffiecient i?^* is defined as 

B NM = ( W + M-j)\N\M\ \ 1/2 mi f, M) (-1) M - ( , 

V 2* +M J m=m ^._ JV) (i-m)!(iV + m-i)!(M-m)!m! 1 j 
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The sums over states above are restricted to Landau levels lying within the pairing width 
around the chemical potential. Using the standard method of evaluating Matsubara sums 
by contour integration we obtain 



/(rii,n 2 ,n 3 ,n 4 ) 



( e -«»n + i)(£ nU + £ n2T )(-£ nii + + &mt) 

M ^ + l)(^ T + £ nii )(£ naT + WX^T - Ut) 

e ^ 4 t + i)(£ + f nii )(f„ 4T - UtXUt + £ 



-i 



+ 

+ 

+ 
-1 



(16) 



The second order term Q 2 which determines the H c2 line agrees, apart from the inclusion 
of the weight function, with the result of MacDonald et a/.il and Rajagopal and Ryan.El The 
six and eighth order terms and fls can also be calculated and they are given in appendix 
[X[ We get the form: 

Q s - n N = ]T aj(T, H)A) + 2 ln-UT, H)A n ■ ■ ■ A H + 

3 h—H 

J2 K h ... j6 (T, H)A n ■ ■ ■ Aj 6 + J2 '/, .,,(•/• H)A n ■ ■ ■ A h (17) 

31—36 31—38 

Thus we have derived the exact quantum mechanical expressions for the expansion coef- 
ficients for Qs — &n up to eighth order assuming a vortex lattice. We have not yet restricted 
the electrons to form pairs with the lowest possible center-of-mass energy (j = 0). The 
result is a multidimensional polynomial in Aj. Going to eighth order permits us to check 
the convergence properties of the series. We could in principle calculate the expansion coef- 
ficients to any order but, as usual, the algebra gets more cumbersome with increasing order, 
and the minimization condition cannot be solved analytically for such high orders. 



C. Self-consistency and minimization of £1$ 

The self-consistent determination of A(r) = V < < 0t( r )V'|( r ) > is equivalent to minimising 
the grand potential with respect to A(r).H In the above formulation, which takes into 
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account the spatial symmetry of the order parameter, this reduces to minimising our multi- 
dimensional polynomial with respect to Aj. Although this is a standard numerical problem it 
is necessary to make further aproximations in order to obtain simple analytical results. The 
instability towards superconductivity is determined by the sign of the expansion coefficients 
aj. Above H c2 we have aj > for all j. The transition to the mixed vortex state occurs 
when one of the aj's becomes negative. The system can then lower its energy by making 
the corresponding Aj nonzero. It has been shown that the instability occurs first in the 
j=0 channel.^ So we have «o < and aj^o > for H ^ H c2 and therefore A 3> Aj^ - We 
can then make the approximation Aj^o = 0, i.e only consider condensation into pairs with 
lowest Landau level center-of-mass motion. We have checked this approximation by solving 
the BdG equation numerically when Aj^ — and when all the Aj's are non-zero. In the 
region of interest there is essentially no difference between the two solutions thus justifying 
our approximation. 

The grand potential now has the Landau form 

aAl + jAt + KAt + nA*, (18) 

(a = ao, 7 = 7oooo etc.) and our self-consistency problem is reduced to a simple one- 
dimensional minimization problem which can be easily solved. To fourth order we have a 
mexican hat potential when we are in the mixed state (a < and 7 > 0) and the minimum 
for the grand potential is obtained for non-zero A Q . Requiring 9a (^s ~ ^v)Ia = gives 

z 3 + o 1 (T, H)z 2 + a 2 (T, H)z + a 3 (T, H) = (19) 

where z = Aq and a\ — a 2 = as = j-. A is the value of A which minimizes 
Qs ~ Equation ( p!9|) is a cubic equation and can be solved exactly. To fourth order we 
obtain 

2^{T,H) 47 (T, H) V ; 

Equation |TU| yields Ao and therefore A(r) and Qs - ^jv as a function of H. The value Ao 
which minimizes Qs ~ will be a function of H and T through the coefficients a, 7, k, rj. 
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Because magnetic quantization has been accounted for exactly, all coefficients and, hence, 
A and Qs ~ &>n are oscillating functions of H for a given temperature T. The condensation 
energy Q$ ~ &>n oscillates 180° out of phase with the normal state Qn close to H C 2. This is 
the origin of the damping of the magnetic oscillations of Qs compared to Qn- The physical 
reason for this effect is rather simple as will be explained in Sec. [IV A[ 

We now consider the magnetization Ms = (c^O^) = (8h(^n + [ps — ^jv]))„- The 
grand potential for a free 2D electron gas Qn can be calculated analytically for the case 
when only two Landau levels are partially occupiedi^S For relatively high T, low H or 
small g-factor, this assumption breaks down but it is then straightforward to calculate Qn 
numerically It should be noted that the chemical potential /x in general is a function of H. 
We have in most of this article, for simplicity, kept the chemical potential \x fixed thereby 
avoiding having to determine \x self-consistently. The oscillatory effect of the chemical po- 
tential is most important for low temperatures (T^0.2) and very clean samples such that 
higher harmonics contribute to the magnetic oscillations. In section [IV (J| we will show that 
even in this case one can to a good approximation consider the chemical potential constant 
in the mixed state. 

III. COMPARISON BETWEEN NUMERICAL DATA AND PERTURBATION 

EXPANSION 

Recently it has been claimed that the degeneracy of the Landau levels should give rise to 
non-perturbative terms in the expression for fls~ &>n making the Gor'kov theory invalid. For 
finite temperature there should be a non-perturbative Ap-term in Eq.fllS|) resulting in many 
interesting thermodynamic effects.il It is therefore of importance to establish the validity of 
the perturbation theory developed in the preceding sections so that we can use it to derive 
results instead of a cumbersome numerical solution. This is essential in the case when many 
Landau levels participate in the pairing since the computation time is very long in this 
regime for the numerical solution. In order to estimate the accuracy of our perturbation 
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expansion, we compare it to an exact numerical solution of the corresponding BdG-equations. 
As mentioned earlier, we have set up a code which solves these equations self-consistently. 
We have chosen parameters such that uje>/uj c = 5, j-^-p = 8.2 and kbT/hto c = 0.28 when 
nf — 12. In Fig. 1 we show the order-parameter Ao as a function of the magnetic field. The 
chemical potential n is fixed. We have plotted both the numerical, the fourth order and 
the eighth order solutions. There is good agreement between the numerical solution and 
our perturbation expansion for both fourth and eighth order. The general behaviour of Ao 
is correctly predicted by both the fourth order and the eighth order expansions. In Fig. 2 
we have plotted the condensation energy Q$ ~ &n- We are measuring energies in units of 
hu c . It is apparent that the contribution Q$ ~ &n has local minima for nf integer. Since 
Qn has local maxima for rif integer the condensation energy oscillates 180° out of phase 
with the contribution from the normal state Qn- We therefore get partial cancellation of 
the normal state oscillations and a damping of the dHvA-oscillations. This is seen in Fig. 
3 where we have plotted the magnetization M = — (9#f2) for both the normal state and 
the mixed state. When the superconducting order starts to increase at nj ~ 10, we get 
significant damping of the dHvA oscillations. Again the agreement with the numerical data 
is good as long as nj £j 12. Eighth order theory tends to agree better with numerical data 
than does the fourth order theory indicating that the perturbation expression is valid. Once 
we go too far into the superconducting state, the perturbation theory starts to disagree with 
the numerical results, also as expected. We see from Fig. 1 and Fig. 2 that the magnitude 
of Ao and Qs — &>n is still fairly well described for n/ > 12, but both fourth and the 
eighth order expansions start to pick up spurious oscillations in the order-parameter and 
in the energy. Q$ ~ &>n actually starts to oscillate in phase with Qn according to the 
perturbation theory. This gives enhancement of the dHvA-oscillations in the mixed state 
as compared to the normal state, as seen from Fig. 3. This is an unphysical effect and is 
absent in the exact solution. Since this enhancement is neither confirmed numerically nor 
experimentally, we conclude that perturbation theory in the single parameter A breaks 
down at this point. It can be shownil that the Gor'kov expansion is convergent if the 
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change in the quasiparticle energies \EQ — (, v \ is not larger than 0(kt,T). We have looked 
at the numerically calculated quasiparticle energies as a function of rif. As expected and 
in agreement with Norman et a/.0 we observe that the quasiparticle bands go from being 
essentially broadened Landau levels close to the transition point to loosing all their Landau 
level structure deeper into the mixed state. For the above specific case we have found that 
for rif ^ 12 the quasiparticle energies are changed so much that the above condition for 
the validity of the Gor'kov series does not hold in large regions of k-space thus explaining 
the breakdown of perturbation theory. We have compared the numerical solution and the 
perturbation expansion for a number of different parameters. Our conclusion is that both 
fourth and eighth order perturbation theories describe well the superconducting state and 
the corresponding damping of the magnetic oscillations near the transition point. However, 
the perturbation theory eventually breaks down when the quasiparticle levels are changed 
too much, in the sense described above. The convergence range of the Gor'kov expansion is 
determined by the temperature k^T. 

The numerical results show total suppression of the dHvA effect once we are deep enough 
into the mixed state. In Fig. 3 the numerical solution shows that M s loses its dHvA structure 
completely for rif ^ 12. This contradicts the recent predictions of a sign shift of the first 
harmonic of the dHvA oscillations.El This prediction is partly based on the assumption 
that the quasiparticle spectrum can be described by a simple splitting of the Landau levels 
into two levels symmetrically placed on around each Landau level even when the actual 
change in energy is rather large (\E^ — ^ v \/hu c m ±0.22). We have found that the low 
lying quasiparticle levels loose their Landau level structure and describe essentially localized 
bound states when the change in energies is of the above magnitude. This crossover to 
localized states makes the argument leading to the sign change of the first harmonic invalid 
and it leads to the suppression of the magnetic oscillations. @ 



IV. DAMPING OF THE MAGNETIC OSCILLATIONS 



15 



A. Physical interpretation 



To get a physical understanding of the superconducting damping of the magnetic oscil- 
lations, it is helpful to consider the ground-state which gives the dominant contribution to 
the grand potential for low temperatures. By analogy with the case of no magnetic field,iil 
our numerical solution is based on the following canonical transformation: 

7kt = Yl [u'jvkflJVkr + ^jvk^-ki] ( 21 ) 

N 

% = £ [ M "k«JVk| - ^kajv-kt] ( 22 ) 

N 

where u v Nk is the coefficient of 0(r) jvk and is the coefficient of <fi*(r) N-k in the Bogoliubov 
amplitudes u(r) and v(r) for the 77'th solution respectively. The corresponding ground state 
of our mean field Hamiltonian is then 

i« r a >«n%^k i i*> (23) 

77k 

where |\& > is a state with all single particle states with energy less that /i — lod empty and 
all single particle states with energy higher than // + u D occupied . We see that Eq. fl23|) 
gives a coherent superposition of states where the _ k |0 > are either occupied 

or unoccupied. When we have a Landau level at the chemical potential /i (jif =integer) it 
does not cost any kinetic energy to make a superposition of states with either occupied or 
unoccupied pairs formed by electrons in that level. The instability towards superconductivity 
is therefore largest when we have \i — (n+ l/2)?ku. Since the grand potential of the normal 
state is at a maximumil when // = (n+ l/2)hu we have that Q$ and fl^ oscillate 180° 
out of phase. This analysis is true for both constant chemical potential and constant number 
of particles. In the latter case one works with the Helmholtz free energy but the conclusions 
are the same. Mathematically the maximum in the damping comes from the fact that when 
the chemical potential is at a Landau level the sum in equation ([LTD is dominated by the 
terms with zero denominators, as an application of 1' Hopital's rule on these terms confirms. 
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Hence a(H) has a local minimum and the superconducting order a local maximum. This 
is the physical picture of the damping of the magnetic oscillations that naturally emerges 
from our formalism. 

Norman et <M interpret the damping of the magnetic oscillations as an effect of the 
broadening of the Landau levels due to superconducting order. An alternate explanation 
has been put forward P. Miller and B. L. Gyorffy.0 that emphasizes the role of non-diagonal 
pairing. There is in fact an intimate link between the two approaches that we now elucidate 
by the following simple calculation: We estimate Q$ — (for simplicity we consider T = 0) 
for the two cases when (I) the chemical potential is at a Landau level (rif integer; maximum 
of the free energy) and (II) when it is exactly between two LL (jif is half an odd integer; 
minimum of the free energy). In both cases we diagonalize the BdG equations approximately, 
but insist on using degenerate perturbation theory, because the diagonal approximation 
breaks down. When n/ is an integer, the lowest lying quasi-particle excitations have the 
orbital character of the n^-st LL, and perturbation theory yields for the quasi-particle 
energy E ^ = \F %\. It can easily be seen that the contributions of the other LL to 
the ground state energy cancel pairwise within degenerate perturbation theory (essentially 
because within degenerate perturbation theory level repulsion is symmetric with respect to 
the unperturbed degenerate level). Therefore the reduction in the maximum of the free 
energy for case (I) in the mixed state is 

k 

to lowest order in the pairing self energy. 

A similar calculation for case (II), when n/ is half an odd integer and the free energy 
is a minimum, gives an energy shift which is of higher than linear order in the pairing self 
energy, because degenerate perturbation theory now leads to complete pairwise cancelling 
for all Landau levels to first order in the pairing self-energy). Therefore the minimum of 
the oscillation is reduced by substantially less than the maximum, which shows that the 
damping of the oscillations is a direct consequence of the broadening of the quasi-particle 

17 



levels accompanied by the mixed orbital character of quasi-particle excitations. 

B. Finite Zeeman splitting 

Inclusion of spin in general reduces the magnitude of oscillations of A and Qs~^n- This 
reduction in the amplitude of the oscillations is due to the fact that spin up and spin down 
electrons now have different energy (unless g = Inm^jm* n = 0, 1,2, . . .). We can never 
have the situation whereby pairing occurs without a cost in kinetic energy. The oscillatory 
effect is therefore damped. The mathematical reason for the reduction in oscillations is that 
for finite spin the numerator in Eq.( [11]) never becomes zero. So we expect the magnetic 
oscillations in the mixed state to be reduced due to spin. The question is whether this 
reduction is larger or smaller than the corresponding reduction in the normal state thus 
giving rise to extra damping effects. For temperatures (or impurity concentrations) such 
that only the first harmonic of the dHvA oscillations is important in both the mixed and the 
normal state £ 0.2hu) c ) and within the region of validity of the perturbation expansion 
of Qs — &n the result is that the amplitude of the first harmonic of the dHvA oscillations 
in the mixed state is reduced by a factor cos(7r|^-). This is the same reduction as in the 
normal state and hence the relative damping due to superconductivity is insensitive to spin 
splitting. This result will be proved in section [V[ We have confirmed this result by solving 
the BdG-equations numerically with and without a finite Zeeman splitting. The reduction 
in the amplitude in both the mixed and in the normal state as compared to the amplitude 
with no spin splitting corresponds very well to a cos(7r|^-) factor in the region where the 
mixed state is described well by the perturbation expansion. Deeper into the mixed state 
the numerical results indicate that the effect of spin is suppressed by the superconducting 
order. The reduction in the amplitude of the magnetic oscillations due to a finite Zeeman 
term is less than the cos(7rr^) factor. This is due to the fact that when the superconducting 
order increases, the pairing interaction starts to dominate the Zeeman term and the effect 
of any finite g-factor is suppressed. 
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So we conclude that within the region described well by our perturbative expansion a 
finite Zeeman term does not alter the rate of the damping of the magnetic oscillations due 
to superconductivity. When only the first harmonic is important the effect of the Zeeman 
term is simply a reduction by a factor cos(7r|^-) for the amplitude of the oscillations in both 
the mixed and in the normal state. Deeper into the mixed state the superconducting order 
starts to suppress the effect of the spin splitting and the magnetic oscillations is less affected 
by a finite Zeeman term. Hence in this region the relative size of the magnetic oscillations 
in the mixed state as compared to the normal state is larger for finite spin splitting and the 
damping is less efficient as compared to the g = case. 

C. Conserved number of particles 

For two-dimensional systems with a fixed number of particles it is well knownH that the 
magnetic field dependence of the chemical potential p>(H) has a strong effect on the magnetic 
oscillations in a normal metal when higher harmonics are important. For low temperatures 
and clean samples the shape of the oscillations look qualitatively different when the chemical 
potential is fixed as compared to when the number of particles is fixed. We have up till now 
mainly considered the case of a constant chemical potential. When the number of particles 
is held fixed we need to consider Helmholz free energy F = Q + Nfi. The chemical potential 
is determined by the equation 

< N >= £ [|< k | 2 /l + K k \ 2 (l - /-«*)] = N (24) 

aNkr/ 

where f\ = (exp(/5i?^ k ) + This is a numerically cumbersome problem since we need 
to solve the BdG- equations self-consistently for a given chemical potential, then calculate 
< N > and repeat the calculation for a new value of fi until Eq . (f2"4|) is obeyed. However it 
is essential that we determine the chemical potential self-consistently. If we naively assume 
that the chemical potential oscillates as in the normal state we would obtain persistent 
magnetic oscillations of the free energy even when the Landau level structure is completely 
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destroyed by superconducting order. In Fig. 4 we have plotted the magnetization when 
the chemical potential is constant (□) and when the number of particles is constant (*) 
for a very low temperature. We have chosen parameters such that cud/^c = 5, j-^-p = 9.0 
and kbT/TiUc = 0.05 and grrf/rriQ = 1 when n/ = 12. For comparison the solid and 
dotted lines give the magnetization in the normal state for rif £ 8.2 for conserved [i and N 
respectively. We see that there is only a significant difference between the two curves close 
to H C 2 (rif ~ 7.7 at H^) when the chemical potential behaves differently in the two cases. 
Deeper into the mixed state the oscillatory behaviour of the chemical potential is damped by 
the superconducting order and it becomes practically constant. This is illustrated in Fig. 5 
where we have plotted rif — ii(H)/Tiuj c — 0.5 as a function of the magnetic field (H c2 ~ 1.5) 
when the number of particles N is constant (solid line) and when the chemical potential 
is constant (dashed line). We see that the oscillations in the chemical potential when N 
is constant are damped in the mixed state. Once the superconducting order has damped 
the oscillations in the magnetization it has also damped the oscillation in fi(H) and the 
behaviour for fixed N is essentially the same as for fixed /i. Thus the conclusion is that 
although there is some difference in the dHvA signal close to H C 2 when N is fixed conserved 
as opposed to fixed \i the overall rate of damping of the oscillations is the same in the two 
cases. 

V. SIMPLIFIED FORM FOR THE DAMPING 

A. The first harmonic of the condensation energy 

To obtain a simple form for the damping, we must take a closer look at the coefficients 
a(H) and 'j(H) given in Eq. (|Tl~D and Eq. (|12"D . As mentioned already, the transition to the 
mixed state occurs when a(H) changes sign. The Gor'kov expansion is most relevant for 
temperatures such that only the lowest harmonics of the dHvA signal are significant. This 
allows us to focus only on the zeroth and first harmonics of the relevant quantities. Thus 
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we take a(H) to have the form: 

a(H) ~ ai(l - H c2 /H) + a 2 cos{2-k^/Tiw c ) (25) 

where ai > and a 2 > 0. The coefficients ai and a 2 will in general depend weakly on the 
magnetic field but we assume they are constant. This is reasonable since for fj,/hu c ^> 1 
the rate of change of a\ and a 2 is very slow as compared to the frequency [imc/Tie of the 
oscillations. The essential physics comes from the sign change of a(H) and its oscillatory 
behaviour, combined with the features of j(H) described below. For simplicity we confine 
ourselves to fourth-order perturbation theory. The fourth-order coefficient j(H) is has the 
form: 

j(H) ~ gi -g 2 cos(2vr/i/^ c ) (26) 

where g\ > and g 2 > 0. Again both g\ and g 2 depend on the magnetic field but this 
dependence is weak as compared to the strong oscillatory behaviour coming from the Landau 
level structure. Note that we have opposite signs for the first harmonics of a(H) and 'y(H). 
In section |V B| we will extract estimates of a 2 and g 2 from Eq. ( |TT| ) an d Eq. ([HJ) whereas the 



approximate expressions for a% and g\ will be given in appendix [BJ Using these approximate 

forms for a(H) and j(H) we get for the condensation energy 

no a 2 (T,H) ( fll (l - H c2 /H) + a 2 cos(27r/i/^ c )) 2 

S N 4r/(T,H) ~ 4( gi -g 2 cos(2irfx/hu c )) ' [ } 

Assuming that g 2 g\ we get the following approximate form for the first harmonic of 

Qs ~ &>n to first order in g 2 /gi- 

2aia 2 92a\ ( , , 2 3g 2 aj 



.(H c2 /H-l)-^(H c2 /H-l) 



cos(27r/i/^a; c ) (28) 



9i 9i ±gi 

where Q(H) n is the n'th harmonic of Q(H). It should be recalled that the above expression 
is only valid for a(H) < 0. When we are deep enough into the superconducting state so 
that we are away from the reentrance region we have a\{H c2 /H — 1) > a 2 . This means 
Za ?l 2 <C — < ^^-(Htf/H — 1) and we can neglect the small constant term 3a2 / 2 . We thus 
get the following form for the first harmonic of the grand potential: 
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(29) 



wherec 



eL x L y H hui c 

2ltTlC 7T 2 



2 A; 6 T _ 2w 2^t 
Arc - — e 



cos(27r/x//ta; c ) 



(30) 



We have written the reduction due to finite temperature in square brackets. 



B. Calculation of the oscillatory terms 



In this section we will derive some approximate expressions for the coefficients 02 and g?,- 
We are interested in how a 2 and g 2 depend on the parameters rif, u D , and T. It turns out 
that it is fairly straightforward to extract this dependence for the oscillatory terms. First 
we note the following approximate identity coming from the law of large numbers: 

gn u n 2 ^ 1 e -("i-n2) 2 /8n 1 ^ 1 c -(ni-n 2 ) 2 /8n f 

where we have assumed that \n\ — n 2 \/ni <C 1 and n\ ~ n/ (i.e u)£>/lo c = 2cr nj). Using 
this formula, Eq. [11], and the Poisson identity we get the following integrals for a(H): 

hu c £ ^^2^(^/2) +tanh(^ 2 /2) w2(ni)w2(n2) = 

ni,n,2 Sni ~l~ £"2 

y e 2, inf (m-i) f dx f dye 2,i(m X -i y) e 4nf tagh(^/2) + tanh(/^/2) w 2 (x)u;2(y) = 

]T/,, m e 2 ™/( m -') (32) 

where ^ = hu c x and w(x) = e -(^/ - 5huJ D) 2 _ e -^ 2 /o- 2 ^ rp Q es ^j ma ^ e j l m w here (l,m) ^ (0,0) 
we write the integral in the form: 

2kf,T f f e ~( x -y) 2 / 4n f e -2(x 2 +y 2 )/a 2 +2ni(mx-hj) 
hm— T ^2 dx ; ; — ; r (33) 

J J ' ^fimj (x - iu' v )(y + 
where io' v = uj v /u c . The first harmonic of a(H) comes from the terms with \l — m\ = 1. 
Taking m = 1 and I = yields the integral: 

r 2Tzix-2x 2 /a 2 f -2y 2 /cr 2 -(x-y) 2 /4n f 

/ dx- — / dy- — (34) 

J x — iu u J y + iu' v 
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We approximate this integral by: 

e 2wix . e -(y-x) 2 /4n f 

dx — / dy — — — (35) 

x — uo' v J y + iu3' v 

since we have assumed 8nf <C a 2 . The integral can be solved exactly and we get : 
J 0il = % 4kb ™ E e- W " ([1 " *Wy/nf)]** /n ' + e <2 [(l - S(f ))>-*"■"' 



u>0 



~ ■ 7r e (36) 

where $(x) = 27r~ 1 / 2 J x e~ s2 <is is the error function and t = (u/ + 27rn/)/^/n/. Here we 
have used that exp(— 27ru; ! ,^o) <S exp(— 27rcov =0 ) for 2n 2 kbT /Tiuj c £ 1, exp(— Am 2 rij) <C 1, 

— 1/2 

and co>„ =0 "C 1. So, in this temperature range the dominant contribution to the first 
harmonic comes from the lowest Matsubara frequency, which makes our approximation 
above self-consistent. The contribution to the first harmonic from the \l — m\ — 1 term given 
l,m ^ can be calculated in the same way; it is proportional to exp(—27r 2 mkT/huj c ) and 
therefore negligible for 2n 2 kbT/huj c ^ J l in agreement with the results obtained by Gruenberg 
et a^After some algebra, the calculations outlined above combined with Eq. (|ll|) lead to 
the following approximate result: 

V2 (t) ^ hT 

The above result that a 2 is proportional to and /c^Te 77 fi "c and independent on 

c<jd is still correct even when a 2 <C 8nj, as long as min(<7, ^Jnj) ^> 1 and 2 r n 2 kbT ' /Tiuj c ~ 1. 
Inclusion of spin is equivalent to making the substitution x — > x + and y — > y — fi^ in the 
integrals Ii >m . This results in a reduction factor cos(7r|^-) in Eq.( |37|) if min( v /ny, a) 

The calculations for gi are very similar to the ones above. Using Eq. ([12]) and the Poisson 
formula we end up with the following integrals determining the dependence of 7 on n/, T 
and c<j£): 

urn . J(xi-X4) 2 + (x3-X2) 2 + {xi-X2) 2 +(x3-X4) 2 ]/8nr 

V ^e^/ft+h-b-U) / dxv . . d ^ x 

££ iU n/ 7 (^-xi)(^ + x 2 )(^-x 3 )(^ + x 4 ) 

v 

~2(x 2 +x 2 .+x'1+x\)/ct 2 2m(l 1 x 1 +lzx i -l 2 X2-hx i )'z:x 1 +X4,, x 2 +x 3 /qoN 
e e "xi+X2,3;3+a;4 v do / 
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where S^'^ is given in Eq. (|i~4"|). Contributions to the first harmonic g 2 come from the terms 
with \l\ + Z3 — 12 — /4I = 1. As in the case for a 2 we can neglect the terms with more than one 
li different from zero when 2ir 2 k\ 1 T /Tiuj c ~ 1. Although we do not have any simple expression 
for we can still extract the dependence on T, u D , and rif. This is because the 

integral over x 2 ■ ■ ■ £4 does not vary appreciably with x\ on a scale ~ lu' u=0 . Using the result 
/ dx ex ^Tx^ f ( x i ) c*- /(0)e~ 27ra/ (a;' > 0) for any well-behaved function f(x) which varies 
slowly for x £ u' and taking l\ = 1, / 2 = ^3 = U = we get the integral: 

khT f e 2mxi ~ e [xl+(x 3 -x 2 ) 2 +x 2 . + (x 3 -x 4 ) 2 ]/8nf 



dx\ — / . . . dx. 



Tlf J %U)' V — X\J '' [vj)' v + x 2 ){iu' v — xz)(iu' v + X4 

: X4,X 2 +X 
-'X2,X 3 +X i 



2(x 2 2 +x 2 +x 2 4 )/a 2 ' = :x 4 ,x 2 +x 3 



The factors 1 / (iu' ± Xj ) in the integrand makes the integral largely independent of any 
long range behaviour determined by a and rif as long as \u'\ <C min(a, s/nj). We therefore 
conclude that g 2 is independent of ujd and that it only depends on rif through the nj 1 - 
factor coming from the four Bq' m coefficients. We also obtain that g 2 is proportional to 
fcfeTexp(— 27r 2 |^). The proportionality constant is found through an exact evaluation of 7 
given in Eq. QI2]). We obtain: 

Again the effect of spin (ie. non-zero g-factor) provide an additional cos(7r|^-) in Eq.( fjQ| ). 
It is not surprising that the oscillatory terms a 2 and (72 are independent of the pairing width 
lod since the oscillations are a consequence of the individual Landau levels going through 
the chemical potential. Likewise the and 1/n/ dependence reflect the fact that the 

probability for two electrons, each with energy (n + l/2)hu c , to form a pair with minimum 
center-of-mass energy is proportional to for high quantum numbers, as can be seen 

from Eq.fl3l"l). This proportionality can be explained via simple phase-space considerations. 
We have tested the dependence of a 2 and g 2 on the different parameters rif, ujb and T and 
we find excellent agreement with our approximate forms. 
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To facilitate comparison with earlier papers we will now formally treat the order param- 
eter A(r) as a free parameter and assume that the oscillatory behaviour of Eq.fllBD only 
comes from the harmonics of the expansion coefficients a, 7 etc. This is of course incorrect 
since the self-consistent order parameter itself is a oscillatory function of the field, making 
the results where the corrections to the harmonics of the dHvA oscillations due to super- 
conductivity are expressed as a power series in of limited validity. However, to 
compare with the earlier predictions we ignore for the moment the oscillations in Ao and 
treat it formally as a free paramter (i.e. (Qs — &n)i = g^Aq — g 2 A^ + ...). Here we focus 
on the A 4 -term since there are discrepancies between the predictions of different authors for 
this term. Since 

A 2 ss (L x L y yiJ dr|A(r)| 2 = ^jQ^K (41) 
we obtain using Eq.(f40|) and Eq.(pO|) the formal result for the fourth order term: 

(n* - n*)il A *-*n» = » n N1 ^- (A) (42) 

StephenB obtained ~ lQVL^ 1 /nf{A/%uj^) A for the same quantity using a different semiclas- 
sical approach. The n/ dependence of the two result agree but the numerical prefactors 
are somewhat different. The above arguments for the n/ dependence of #2 can easily be 
generalised yielding that the nf dependence of the first harmonic of the A 2n -term is nf~ n l 2 . 
This rif dependence agrees with the result obtained by Stephen whereas it disagrees with 
the 3//2 -dependence for the A 4 -term obtained by Maniv et a/.H We cannot overemphasize 
the fact that the above scheme to calculate the damping of the oscillations due to super- 
conductivity is incorrect, since it ignores the oscillations in A as a function of the field. To 
include those we have to use a self-consistent order parameter and hence Eq. fl28|) . 
One debated issue is the possibility of reentrance for type-II superconductors. 

The 

oscillatory behaviour of a due to the Landau level structure gives rise to the possibility 
of several solutions to a(H) = for a given temperature. This should be reflected in a 
highly oscillatory behaviour of the transition line H C 2(T,H). Such an oscillatory behaviour 
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has never been observed experimentally. Using the approximate expressions for a\ and ci2 
we can estimate the temperature below which there is reentrance and such oscillations in 
H C 2 should occur in a two dimensional metal. We obtain that when there is no impurity 
scattering, no Zeeman splitting, and rif ~ O(10 2 ) one should observe these oscillatory effects 
in H c2 in a 2D metal for temperatures lower than kbT/hu) c ~ 0.3. However, inclusion of spin 
reduces the amplitude of the oscillations of a by a factor cos(7r|^-) close to the transition 
line. Assuming that impurities reduces the oscillations by a factor exp(—27i 2 kbTr)/huj c ) 
where To is the Dingle temperature we obtain that there will not be any reentrance if 
kbTo/fiuj c ~ 0.2 no matter how low the temperature is. In the case of the experiments being 
done on k — (ET) 2 Cu(NCS)2^ the experimental parameters are such that k^T^/huc « 0.27 
and |cos(7r|^-)| ~ 0.3. They will therefore never observe these reentrance effects. The 
magnetic oscillations in the thermodynamic quantaties will of course still be there since a 
and 7 are still oscillatory. 



C. Approximate results for damping 

In this section we will draw some conclusions from the general form of the damping of the 
dHvA-oscillations due to the growth of the superconducting order described by Eq. (|28l). The 
first thing we notice is that in this approximation the superconducting damping has a simple 
polynomial form in (H C 2/H — 1). The damping is maximum for {H&jH — 1) = For 
(H C 2/H — 1) > 2^ the damping decreases when we go deeper into the superconducting state 
and for (H c2 /H - 1) > ^£1 the magnetic oscillations are enhanced by the superconducting 
order. This explains the observations made in section |TJ. The in-phase oscillations between 
the fourth order Qs — and Qn are due to the oscillatory behaviour of ^(H). Since ■j(H) 
oscillates in phase with Qn we will get the enhancement of the oscillations of Q,s compared 
to fijv when the smooth part of a(H) is sufficiently large. Again we must emphasize that 
this is obviously a sign that our perturbative scheme has broken down and does not reflect 
any physical effect. 
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To make any quantitative predictions we need to use our approximate expressions for 
and <?j. Since we only have very good approximations for a 2 and g 2 and for the tem- 
perature and spin dependence of a\ and g\ we will concentrate on properties that can be 
derived from these results. From Eq. fl28|) and the temperature dependence of dj and gi 
we conclude that the first harmonic of the condensation energy (Q$ ~ &n)i is proportional 
to A^Texp(— 27r 2 |^). Since we also have &ni oc fc&Texp(— 27r 2 |^) this means that the 
magnetic oscillations have the same temperature dependence in the mixed state as in the 
normal state. This result agrees with the general theory (see Schoenbergil Sec. 2.5 and 
Sec. 2.3) valid for any part of the grand potential which is proportional to cos(n/hu c ). It 
is also confirmed by experimental observations.0 Likewise the effect of spin on (Q$ ~ &n)i 
is a reduction in the amplitude by a factor cos(7r|^-). This is the same reduction factor as 
for the oscillations in the normal state.il We thus have no extra damping effects due to spin 
close to the transition line where the perturbation theory is valid. 

We can now examine whether the arguments based on the Gor'kov expansion leading to 
a sign change of the first harmonic are valid. Naively one would expect a sign change since 
the contribution from the condensation energy to the magnetic oscillations is in antiphase 
with the normal state oscillations. When the system is deep enough into the mixed state 
the superconducting oscillations would dominate leading to a sign change of the magnetic 
oscillations. Extrapolating the rate of the damping close to H c2 obtained from the Gor'kov 
expansion Maniv et a/El have estimated the magnetic field H inv < H c2 at which this sign 
change should occur. We are now able to show that this argument based on the perturbative 
expansion of the grand potential is incorrect. From Eq. fl28|) we obtain that the maximum 

a 2 

amplitude of the antiphase oscillations of £ls — is given by j^. Using our approximate 
expressions for a 2 and g 2 we get 

a\ L x L y a 2 x n - 2 ^h? . , 

— * -^-hTe (43) 

Comparing this amplitude with the contribution from the normal state oscillations given in 
Eq. (|30|) we see that our perturbation scheme roughly predicts a maximum damping of 50%. 
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It must be emphasized that this does not mean that the damping of the model described 
by the Hamiltonian in Eq.( |3]) has a maximum of 50%. However, using the result above 
combined with the results in section |I1 we can conclude that neither the argument based 
on the Gor'kov expansion nor the arguments based on a simplified form for the quasiparticle 
spectrum leading to an inversion of the first harmonic of the dHvA signal are valid. 

VI. COMPARISON WITH EXPERIMENT 

In this section we present a typical result for the damping of the magnetic oscillations 
obtained from our theory when many Landau levels participate in the pairing. We have 
chosen parameters such that we can compare our result with the experimental observations 
made by van der Wei et aM First we compare our approximate expressions for the damping 
from Eq.(p8|) with the result based on the exact evaluation of a and 7 from Eq. ([□]) and 
Eq. (|I2"D. We used the a set of parameters such that k b T/hu c = 0.25, = 2.315, and 
uj£,/uj c = 75 when rif = 175. There is no Zeeman effect and the chemical potential is 
conserved. In Fig. 6 we have plotted the magnetization for both the normal state and the 
mixed state calculated from the perturbative expansion to fourth order as a function of 
rif. The perturbation theory predicts a substantial damping of the oscillations over many 
periods reaching a maximum for rif ~ 170. At the maximum the first harmonic is damped 
approximately 50 % in agreement with the result in the previous section. As we go deeper 
into the mixed state, the damping decreases according to the perturbative scheme. Based 
on the results in Section |TJ, we expect the perturbation theory to describe the damping 
well for n/^.170 . Due to the large number of Landau levels involved in pairing, we have 
not undertaken the exact numerical calculation for this set of parameters. In Fig. 7 we have 
plotted M s calculated from the exact evaluation of a(H) and j(H) and calculated from 
Eq. (p8|) . We see that the simplified expression reproduces the perturbative predictions well. 

The above parameters approximate the experiment performed by van der Wei et a/.i 
on the essentially 2D organic superconductor k— (ET) 2 Cu(NC 'S) 2 - To compare with the 
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experimental data we will formulate our results in terms of a field dependent quasiparticle 
scattering rate r defined such that e~ 7T ^ ulcT gives the damping of the first harmonic of the 
dHvA oscillations due to superconductivity. From Eq. we get: 

r -i = _^ in(i + (n 5 - n^i/n^) ~ - UeP ^ (h c2 /h - 1) (44) 

where we have used Eq. (p8|) . The approximate equality is only valid for -C i? C 2 /H — 1. 
Using the expressions for Oj, and we can now compare this expression with the 
experimental observations. Unfortunately the experimental value of H c2 is uncertain. The 
transition from the normal state to the superconducting state occurs over a field range of 
approximately 2T.0 This gives a 'smooth' variation of the r _1 on entering the mixed state 
which our theory cannot account for. To model this transition region we use the method 
introduced in ref. 2 by including a Gaussian spread in H c2 In Fig. 8 we have plotted the ex- 
perimental data for r -1 (bars) measured in (THz) as a function of 1/B measured in Tesla -1 . 
The solid line is our theoretical prediction based on Eq. ( [44]) including a Gaussian spread in 
H c2 . The agreement between theory and experiment is good. It should be noted that we 
have no fitting parameters apart from H c2 . However, without a more reliable measurement 
of H c2 a precise comparison between our theory and the experimental observations cannot 
be made. 

VII. CONCLUSION 

In this paper we have examined the dHvA oscillations in the mixed state of a type II 
superconductor in the 2D limit using both a numerical solution of the BdG equations and an 
analytical theory based on a self-consistent Gor'kov expansion. The use of translational and 
rotational symmetry has simplified the analysis such that we have been able to calculate the 
expansion coefficients exactly to any order without using semiclassical or other approxima- 
tions. Comparison with the exact numerical solution has showed that perturbation theory 
works well close to H c2 thereby disproving recent claims of non-perturbative effects. We 
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have found that the condensation energy oscillates in antiphase with the normal grand po- 
tential, thus producing damping of the dHvA oscillations in agreement with numerical and 
experimental results. The damping is directly connected with the enhancement of super- 
conductivity when we have a Landau level at the chemical potential. We have excluded 
the possibility of a sign change of the first harmonic of the dHvA oscillations in the mixed 
state. The effect of spin and a conserved number of particles as opposed to a conserved 
chemical potential was examined. Using a simple approximate form of our analytical theory 
valid when many Landau levels participate in pairing we have compared our theory with an 
experiment on the quasi 2D organic superconductor k— (ET) 2 Cu(NC 'S) 2 . We have found 
good agreement. However, due to experimental uncertainty about H c2 any quantitative 
comparison is impossible. 
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Using the symmetry of the vortex lattice and making the restriction A^ = we obtain: 
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(Al) 



where 



f(n 1 ,n 2 , . . . ,n 2 z) 



1 



-1 



(A2) 



and 
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- ^ Xi 1 (k)---X i! (k)x* +1 (k)---X* 2i (k) (A3) 



Xj(k) = v 7 / ]T e^^e-^Vi^y + fc^//)) (A4) 



Likewise the eighth order term gives for Aj^o = 0: 



n n -°0 -°0 ^0 -°0 "n 7 +n 6 ,n5+n4,n 3 +n 2 ,ni+n g 

The Matsubara sums and the k-sums can be calculated as in the fourth order case. 

APPENDIX B: 

In this appendix we will extract the dependence of a\ and g\ on nj, T, a and spin. 
This is considerably harder than for a<i and g<i because we do not have any oscillatory 
factor in the integrals that would make the long range behaviour of the remaining integrand 
insignificant. It turns out that it is still fairly straightforward to derive the temperature and 
spin dependence of a\ and gi, whereas we have to make some rather drastic approximations 
to obtain the dependence on n/ and a for g\. 

The smooth part (zero harmonic) of a(H) comes from the terms Iij in Eq.(^). We first 
look at the term I = m = 0. Making the variable substitution v = u = ^7= we get the 
following integral: 

J , = ° I du ( dve~ ( ^ +2)u2 \t&nh[Ka{v + u)} + ttmh[Ka(v - it)])- (Bl) 

y/ZTTUf J J V 

where K = /3huj c /2\f2 S> 1 determines the temperature dependence of the integral. Since 
K is only important around the region v ~ which does not contribute significantly to the 
integral, we conclude that J ,o is independent of the temperature to a very good approxima- 
tion. Since similar calculations to the ones in Sec. |V B| show that for 27r 2 kbT/hu! c ^ J 1 we can 



neglect the contribution to the zero harmonic from the J/^-terms where I 7^ 0, we conclude 
that a\ for is independent of the temperature for temperatures that are not too low. We 
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have checked this independence against the exact result given in Eq.(|TT|) and found very 
good agreement. To obtain the dependence on nj and a we make the simplification 

if Id < Id 

temh[Ka(v + u)} + tanYi[Ka(v - u))v~ l ~ < (B2) 

r, if \v\ > \u\ 
{ \v\ ii ii 

which is a very good approximation since K ^> 1. It is exact for T = 0. The integral can 
be solved and we obtain: 



4 

^o,o ^ , In 



a 2 



+ 1 + 



o 2 \ ,a 2 



n f /a 2 ) n f 

where we have assumed a 2 ^> Arif. This yields the result 



+ 2 ~ 21n(— ) (B3) 



Y 2 \^ ax 

1 4nL x Lyl 2 hu c ^ 

The expression for ai is independent of any spin effects for min(y / n/, a) ^> g. We have again 
checked the independence of on rif, a, and spin against the exact result and we find very 
good agreement. 

The dependence of g± on nf, a, and T is determined by the integrals in Eq.(^) for 
which h — h + h ~ h = 0. Again it turns out, that for 27r 2 /cfoT ' /Tiuj c ^ 1 we can neglect the 
contribution to g\ from the terms with l\ — l 2 + h — U = and max(|/i|, |Z 2 |, |^|) > 0. 
Using Eq.(|l^) we can rewrite the integral with / x = l 2 = l 3 = U = as 



J dx\. . . dx^ 



twdh.{Kxi/2) tanh(-FCz 3 /2) 



£3 — ^1 V ( X l + X 2)(%1 + X4) (x 3 + X 2 )(x 3 + X4) i 

1 / tanh(iTx 2 /2) tajah(Kx 4 /2) 



+ 



xa - x 2 \(xi + x 2 )(x 2 + x 3 ) (s 4 + Xi)(x4 + x 3 

- J [{xi-X4) 2 +(x3-X2) 2 +(x 1 -x 2 ) 2 +{x3-X4) 2 ]/8n f -2(xl+x^+xl+xl)/a 2 ^x 1 +x 4 ,x 2 +x 3 



X 



^Xl+X 2 ,X 3 +X4 



(B5) 



where K = huj c /kbT determines the temperature dependence. Again for min(y / ny, a) ^> g, 
gi will be independent of spin effects. As in the case of ai, it is fairly straightforward to 
see that since 1/K -C min(y / n/, a), the integral and therefore g\ are independent of the 
temperature to a very good approximation. We have checked this independence against the 
exact result given in Eq . (|T^) and find very good agreement. 

32 



To make any progress in determining the dependence of g\ on a and rif we need some 
simple expression for S^^*'^^. As a rough approximation we make the following simpli- 
fication to: 

■^n 1 +n 4 ,n 2 +n 3 / c c \ L x L y 

"ni+n2,n3+7i4 \ u ni,n-$ ' U n2jrM J v / 

This is based on the fact that x Jl (k)x* 2 (k) in general is a complex number for ji ^ j 2 . 
When the k-sum is performed the phase factor will change 'randomly' and make the sum 
approximately zero. Physically it corresponds to ignoring cases where electrons in four 
different Landau levels interact. Using this simplification, Eq. (^Tj) and the Poisson formula 
we get from Eq. [1^ the following integral determining the dependence of gi on rif and a: 

where uo 1 = uj v /uo c . We have again assumed 2ir 2 k} ) T/hLJ c ~ 1. Assume now that 8rif <C a 2 . 
We approximate the integrals by: 

L> ^ / g / dx 2 —— ] / dx 4 — — (B8) 

rif J (ioj — x) J iu' + x 2 J too' + x 4 

We will now show, that in this approximation the sum of the integrals is largely independent 
of rif and therefore g\ oc 1/nf. The integral can be solved and we get: 

hT f 



= ^ [1 " W + ^ - * ( v# )|2e 2 "' (B9) 



'7 

where we again have = -7= Jq dte~ f2 . Eq. (|B7|) can then be written on the form 



2nnf 



Ax^ (-™ n e x * l2 [l - H^)} + V2^j vr 2 [l - $(^)]V"/ 2 



(BIO) 



where x n = an d Ax = %^ k ^j - Since Ax < 1 we can approximate this sum by an 

integral and we therefore conclude that g\ is independent of the temperature in agreement 
with the result above. Furthermore we obtain g\ oc I /rif for rif large and g\ independent 
of a. When a 2 8rif does not hold the calculation is the same as above. We just have to 
substitute l/4n/ with l/4n^ + 2/ a 2 in the integrals. The 1/rif dependence coming from the 
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Bq 1 n factors in Eq. [12] is unaltered and we still get that g\ oc 1/rif for min(y/ny, a) large 
and that g\ is independent of a and the temperature. By calibrating g\ through an exact 
evaluation based on Eq. (|i~2"D , we obtain: 

91 ~ {L x L y fP{hu c fbAn f (BU) 

where gi is defined in Section [V Aj It should be noted that the dependence of g\ on nj 
and a in the above expression is only approximate and rests on the various simplifications 
made. We have tested the above expression against the exact result and we find that the 
dependence on nf and a fits to an accuracy of 20%. 
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Figure Captions 



Fig. 1: The order parameter A vs rif calculated numerically (solid line), to fourth order in 
Ao (dashed line), and to eight order in Ao (dash-dot line). 

Fig. 2: The difference Q,s ~ &>n in the grand potential between the mixed state and the 
normal state. The solid line is a numerical calculation, the dashed line is fourth order 
perturbation theory, and the dash-dot line is eighth order perturbation theory. 

Fig. 3: The magnetization vs n/. The solid line is a numerical calculation, the dashed line 
fourth order, the dash-dot line eight order, and the dotted line is the magnetization in the 
underlying normal state. 

Fig. 4: The magnetization when the chemical potential is constant (□) and when the number 
of particles is constant (*) for a very low temperature. The solid and dashed lines are 
the normal state magnetization for fixed chemical potential and fixed number of particles 
respectively. 

Fig. 5: as a function of the magnetic field for fixed chemical potential (dashed line) and 
fixed number of particles (solid line). 

Fig. 6: The magnetization vs rif in the mixed state (solid line) and in the underlying normal 
state (dashed line). 

Fig. 7: The magnetization vs rif in the mixed state. The solid line is the first harmonic of 
the perturbative calculation and the dashed line is obtained from Eq. fl28|). 
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Fie. function of 1/B. The solid line is the theoretical prediction and the bars 

are the experimental data.i 
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